Brownian motion of a self-propelled particle 
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Overdamped Brownian motion of a self-propelled particle is studied by solving the Langevin 
equation analytically. On top of translational and rotational diffusion, in the context of the presented 
model, the "active" particle is driven along its internal orientation axis. We calculate the first 
four moments of the probability distribution function for displacements as a function of time for 
a spherical particle with isotropic translational diffusion as well as for an anisotropic ellipsoidal 
particle. In both cases the translational and rotational motion is either unconfined or confined to 
one or two dimensions. A significant non-Gaussian behavior at finite times t is signalled by a non- 
vanishing kurtosis "/{t). To delimit the super-diffusive regime, which occurs at intermediate times, 
two time scales are identified. For certain model situations a characteristic t'^ behavior of the mean 
square displacement is observed. Comparing the dynamics of real and artificial microswimmers like 
bacteria or catalytically driven Janus particles to our analytical expressions reveals whether their 
motion is Brownian or not. 

PACS numbers: 82.70.Dd, 05.40.Jc 



I. INTRODUCTION 



There are numerous realizations of self-propelled par- 
ticles [l|, 0] in nature ranging from bacteria |3l-[l0{ and 



spermatozoa [lll - ll3j to artificial colloidal microswim- 
mers. The latter are cither catalytically driven uM- 
Il9| or navigated by external magnetic fields |20l - |22| . 
but also biomimetic propulsion mechanisms can be ex- 
ploited [231 . On the macroscopic scale, vibrated polar 
granular rods |24l - [26j and even pedestrian s l2'i1 provide 
more examples of "active" particles [2^ [23 • A suit- 
able framework for theoretical modeling of self-propellers 
is provided by the traditional Langevin theory of an 
anisotropic particle with translational and orientational 
diffusion including an effective internal force [3Ct |3l| in 
the overdamped Brownian dynamics [3 21 |33|. The di- 
rection of the theoretically assumed internal propulsion 
force (corresponding to an imposed mean propagation 
speed) fluctuates according to rotational Brownian mo- 
tion |34| - [36J . It is a challenging question whether real 
self-propelled particles can at least in a rough way be 
covered according to this simple Brownian picture. While 
for "passive" ellipsoidal particles a comparison revealed 
very good agreement with the picture of Brownian dy- 
namics [37h39| , this has never been undertaken for self- 
propellers. Any deviations point to the relevance of hy- 
drodynamic interactions, non-Gaussian noise, or fluctu- 
ating internal forces which are beyond simple Brownian 
motion. 

Despite its simplicity, the Brownian motion of 
anisotropic particles pOl - l43 | has only been considered 
in the absence of internal driving forces either in the 
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bulk [371, |43J or in an external force field derivable from 
a potential [4J|. For "passive" rodlike particles the 
Smoluchowski-Perrin equation [SJ, Ea] has been solved 
exactly in two IdGf as well as in three [43] dimensions. 
With regard to self-propellers, so far analytical results 
are only available if the orientation vector is confined to 
two dimensions. For rodlike particles [4^, H^] the first 
two [23], and for spherical particles the first four [50] 
moments of the probability distribution function for dis- 
placements were calculated. In this paper, we close the 
remaining gaps by presenting a comprehensive model and 
calculating the first four moments of the displacement 
distribution function for all relevant situations. First, we 
provide analytical results for an anisotropic self-propelled 
Brownian particle in two dimensions. Furthermore, both 
the situations of an isotropic and of an anisotropic parti- 
cle are extended to the full three-dimensional case where 
the orientation vector is unconfined. 

Studying the mean square displacement reveals a 
super-diffusive regime at intermediate times, which is 
characterized by a t^ time dependence for most cases 
and beyond that by a t^ behavior for some special cases. 
Moreover, two time scales that delimit the super-diffusive 
regime are identified. These can be extracted from the 
results for the mean square displacement or from the nor- 
malized fourth cumulant (kurtosis) j{t) of the probabil- 
ity distribution function for displacements, which mea- 
sures the non-Gaussian behavior as a function of time 
t. For small and very large times, the kurtosis vanishes 
indicating Gaussian behavior, but due to both particle 
anisotropy and self-propulsion, j(t) is non- vanishing for 
intermediate times. While Han and coworkers [37[ found 
the kurtosis of "passive" particles to be positive for i > 
with a simple maximum at finite time |37| . here we find 
that a propulsive force tends to make the kurtosis nega- 
tive. There is a rich structure in ^{t) revealing different 



non-Gaussian behavior at different time scales. Quite 
generally, the propelling force induces a negative massive 
long-time tail in j{t) which tends to zero as 1/t. This 
prediction can in principle be verified in experiments on 
self-propelled particles. 

The paper is organized as follows: In Sec. HI] we present 
and motivate the various model situations that are con- 
sidered in Secs. lIIll to lVII of this paper. In each case the 
first four displacement moments are calculated analyti- 
cally and the results are analysed based on appropriate 
figures. Finally, we conclude and give an outlook on fur- 
ther expansion of our model in Sec. IVIIl 



II. REMARKS ABOUT THE VARIOUS MODEL 
SITUATIONS 

In this section, we give an overview of the situations 
to which the model is applied in this paper (see also 
Fig. [T]). In general, the model consists of an isotropic or 
anisotropic self-propelled particle which undergoes com- 
pletely overdamped Brownian motion. To describe the 
propulsion mechanism on average, we theoretically as- 
sume an effective internal force F = Fu that is included 
in the Langevin equation. The orientation vector u is 
introduced to specify the direction of the self-propulsion. 
Depending on the number of translational degrees of free- 
dom, in some of the cases to be covered this force is 
projected either onto a linear channel or onto a two- 
dimensional plane. To characterize the different situa- 
tions depending on the number of degrees of freedom of 
the particle, we introduce the following notation: The 
(D,^, CT)-model refers to the situation with D transla- 
tional degrees of freedom and d orientational degrees of 
freedom. The possible values for these parameters are 
D G {1,2,3} and d S {1,2}. The parameter a e {s, e} 
refers to the shape of the particle. While a = s relates 
to a spherical particle, for an ellipsoidal particle a = e 
is used. When no specific value is given for one of these 
parameters, we refer to the group of models with an ar- 
bitrary value for that parameter. 

We will first refer to the {D, 1, s)-model, which is de- 
picted in Fig. [TJa) (theoretical investigation in Sec. IIIip. 
This system consists of a self-propelled spherical par- 
ticle whose rotational motion is constrained to a two- 
dimensional plane. To study the behavior of the parti- 
cle, we first refer to the one-dimensional translation in 
x-direction fSecs. HITXl to HlTC]) . In experiment, this sit- 
uation can be achieved by confining swimmers by means 
of external optical fields [51, 52} , for example. After in- 
vestigating this (1, l,s)-model the results can easily be 
transferred to the two-dimensional case ((2, 1, s)-model), 
which is done in Sec. IIIIDI This model situation is es- 
pecially useful to describe the motion of a self-propelled 
particle on a substrate. 

As the assumption of spherical particles is not justifi- 
able in many experimental situations, the model is gener- 
alized to ellipsoidal particles (see Fig.[IJb)) in Sec. lIVI bv 
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FIG. 1. (Color online) Sketch of the various situations to 
which the model is applied: a) (D, 1, s)-niodel, b) (D, l,e)- 
model, c) {D, 2, s)-model and d) {D, 2, e)-model. The nota- 
tion is explained in the text. For D = 1 in subfigures (a) and 
(b) and for _D < 3 in subfigures (c) and (d) , the effective driv- 
ing force along the particle orientation u is projected onto the 
respective number of translational dimensions. To cover the 
(1, 1, s)-model, for example, in subfigure (a) only the motion 
in the x-direction is considered. 



investigating the (Z?, l,e)-model. The more complicated 
coupling between rotational and translational motion as 
opposed to spherical particles leads to qualitatively dif- 
ferent results. 

Besides particles moving on a substrate with one ori- 
entational degree of freedom, we study freely rotating 
self-propelled particles. Here, the case of free transla- 
tional motion in the bulk ((3,2,a)-model) is interesting 
as well as situations in which the translation of the par- 
ticle is constrained either to a linear channel ((l,2,cr)- 
model) or to a two-dimensional plane ((2, 2, cr)-model). 
Again, a spherical particle (see Fig. [TJc)) is discussed 
first (Sec. |V|, before the most general case of a freely 
rotating ellipsoidal particle (see Fig. [IJd)) is considered 
in Sec. IVD 



III. SPHERICAL PARTICLE WITH ONE 
ORIENTATIONAL DEGREE OF FREEDOM 

This section contains the theoretical considerations 
concerning the situation in Fig. [TJa) ((£), 1, s)-model). 
Although a spherical object is regarded here, one has to 
consider that a certain direction is specified through the 
theoretically assumed driving force F = Fii. The two- 
dimensional motion of the self-propelled particle can be 
described by the coordinates x and y of the center of mass 
position vector r(i) = {x,y) and the angle (j) between e^ 



and u — {coa(j),am(j)). Thus, the basic Langevin equa- 
tions are given by 



^ = (3Dt [Fu ~WU + i] 

— - = PDrg ■ e^ . 
at 



(1) 
(2) 



Here, f(t) and g(i) are the Gaussian white noise ran- 
dom force and torque, respectively. They are character- 
ized by (Mt)) = 0, {h{t)f,{t')) - 2%5(i-i')/(/3'A), 
{g.{t)) = 0, and {g^{t)gAt')) - 2<5„^(i - t')/{P^Dr), 
where the indices i and j refer to the respective compo- 
nents, 5ij is the Kronecker delta, and (. . .) denotes a noise 
average. C/(r) is an external potential. The prefactors in 
Eqs. (HI) and ^ consist of the inverse effective thermal 
energy /3 = {k^T) on the one hand and the transla- 
tional and rotational short-time diffusion constants Dt 
and Dr on the other. In the case of a spherical parti- 
cle with radius R these two quantities fulfill the relation 
Dt/Dr = 4i?^/3, which is used in the following analytical 
expressions for the displacement moments. The Langevin 
equation ([2]) can easily be derived from the more general 
vector equation (du)/(dt) = (3D,.g,{t) x u. 

As is a linear combination of Gaussian variables ac- 
cording to Eq. (m , the respective probability distribution 
function has to be Gaussian as well and proves to be 
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where 0o = 0(^ = 0) is the initial angle. For the theoreti- 
cal analysis, the two-dimensional motion in the xy-plane 
can be split up into its components in x- and y-direction, 
respectively. In the following subsections, we first refer 
to the z-component of Eq. ([1]) . As some calculations for 
the (1, 1, s)-model have already been presented in 15 0| in 

I 



more detail, we only summarize the most important re- 
sults and briefly refer to systems with additional linear 
or quadratic potentials after that. 



A. The (l,l,s)-model 

Integrating the averaged Eq. ([T]) for [/ = over time 
and considering only the x-component yields 



{x{t) - xo) = -PFR^ cos(</)o) [1 - e- 



(4) 



and 
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for the mean position and the mean square displacement. 
As usual the skewness S and the kurtosis 7 are defined 
as 



S ■ 



and 
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respectively. A non-Gaussian behavior is manifested in 
non-zero values for S and 7. Using the notation F* — 
(3RF for spherical particles and a scaled time r — Djt, 
the third and fourth moments are given by the analytical 
results 
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FIG. 2. (Color online) Mean square displacement of a spheri- 
cal particle with initial orientation angle (f>o = O.Stt for differ- 
ent values of F* — jSRF. The different time dependences in 
the various regimes are illustrated by straight lines. 




FIG. 3. (Color online) Mean square displacement of a spher- 
ical particle with effective force Fs — 100 for different values 
of (pQ. The regime of super-diffusive motion is defined by the 
time scales ii and t2. While t2 is independent of the initial 
angle <j>o, the time scale fi is the bigger the more the initial 
orientation of the particle deviates from the a;-direction. 



In the case of large forces (3RF ^ 1, the particle mo- 
tion (see Figs. [5] and [3]) is separated into three quali- 
tatively different time regimes, two diffusive regimes at 
short and at large times, and a super-diffusive regime at 
intermediate times, which is characterized by a f^ or a 
t^ behavior of the mean square displacement, depending 
on the initial particle orientation. These regimes are sep- 
arated by the two time scales ti and t2, respectively: At 
early times i <C ii, the particle undergoes simple transla- 
tional Brownian motion, which is governed by the short 
time translational diffusion term {8/3)R^Drt in Eq. ([5]). 
As seen in Figs. [2] and [3] the mean square displacement 
displays a crossover to an intermediate super-diffusive 



regime at a time scale ii , which, in turn, depends on the 
initial orientation 0o and the effective force F* — /3RF. 
In particular. 



t 



{■i/2){PRFDrr\ 
(3/2)[cos(0o)/?i?-F]~ 



'DZ 



cos(0o)| < l/V^RF 
cos(0o)l >1/VPRF. 
(10) 
If the initial orientation has a sufficiently large compo- 
nent parallel to the z-axis, i.e., if |cos(0o)| > l/y//3RF, 
the mean square displacement displays a crossover to a 
ballistic regime, which is governed by a scaling relation 
{{x{t) — xq)'^) cx t^. The crossover is observed the ear- 
lier the larger the inital force component | cos{(j)o)\/3RF 
parallel to the x-axis. On the contrary, if the initial par- 
ticle orientation points along or almost along the y-axis 
at i = 0, i.e., | cos(0o)| < ^/VPRF: the crossover time ii 
is substantially larger. In the latter case ti is the time it 
takes the particle to undergo an angular displacement by 
rotational diffusion, such that the projected force onto 
the X-axis becomes as large as demanded in the former 
case. Only then does the force bring about a lateral 
displacement that is compatible or larger than the dis- 
placements due to original translational Brownian mo- 
tion. Due to this multiplicative coupling of a diffusive 
and a ballistic behavior for the angular and the transla- 
tional displacements, respectively, the mean square dis- 
placement shows a super-ballistic power-law behavior for 
t > ii, with {{x{t) - xo)^) cx t^. 

The intermediate regime is terminated by free rota- 
tional Brownian motion at the second time scale ^2 — 
D^^, beyond which the particle motion is diffusive again. 
As already reported in Ref. [50|, the long-time transla- 
tional diffusion constant is given by 
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(11) 



As can be seen in Fig. lU the beginning of the super- 
diffusive regime also shows up in the kurtosis. Here, the 
deviation from zero clearly indicates the crossover to non- 
Gaussian behavior. Interestingly, the kurtosis features a 
pronounced long-time tail. Therefore, the behavior of the 
particle is still non-Gaussian when its motion is (nearly) 
diffusive again. Analysing the analytical result for the 
kurtosis gives the leading long-time behavior as 
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As the amplitude vanishes for F* — 0, this negative 1/i 
long-time tail in 7(i) proves to be characteristic for self- 
propelled particles. 

The previous calculation of the displacement moments 
can also be done for systems in which the self-propelled 
Brownian particle is exposed to x-dependent linear or 
quadratic potentials. This is illustrated in the following. 
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FIG. 4. (Color online) Kurtosis 7(t) of the probability dis- 
tribution function ^(a;,i) of a spherical particle with initial 
orientation angle 0o = O.Stt for the same values of Fl = PRF 
as used in Fig. [2] Comparing Figs. [2] and |4] shows that the 
time scales tia, tn, and tic can be extracted from the plots of 
the kurtosis as well as from the mean square displacement. 



FIG. 5. (Color online) Mean position of a spherical parti- 
cle with driving force F^ = 10 that is exposed to an exter- 
nal square potential. The strength of the potential U{x) = 
{l/2)kx^ is determined by the parameter A = {A/i)l3kE? . 
The dimensionless quantity c = xo/R is the distance between 
the initial position of the particle and the position of minimal 
potential given in units of particle radius R. 



B. The (1,1, s)-model with an additional linear 
potential 



The Langevin equations in this case are obtained from 
Eqs. H]) and ^ by simply inserting a linear potential of 
the form U{x) ~ mgx into the first eomponent of Eq. ([l}. 
Thus, the motion of a particle that is exposed to gravity 
is described by the moments 
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While the first moment (Eq. ([T^ ) is a simple superposi- 
tion of the terms due to the self-propulsion of the parti- 
cle and the external force, respectively, the mean square 
displacement (Eq. (ITil) ) has an additional term which de- 
pends on both of these forces. 



FIG. 6. (Color online) Mean square displacement of a spher- 
ical particle that is exposed to an external square potential. 
The plots refer to the same parameters as the plots in Fig. [5] 



C. The (1,1, s)-niodel with an additional square 
potential 

Based on the Langevin equations ([1]) and ([2]) one also 
obtains analytical results for the harmonic oscillator or 
square potential U{x) = {l/2)kx'^. Using the dimen- 
sionless parameter A = {4:/3)l3kR'^, vifhich determines the 
strength of the square potential, the mean position is 
given by 
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and the mean square displacement is calculated as 
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Figure [5] shows that the mean position of the particle 
reaches the position of the minimal potential and stays 
there. Due to the square potential, the mean square dis- 
placement (see Fig. ^ does not diverge as in the cases 
that have been regarded up to this point. 



D. The (2, l,s)-model 

In this section we briefly want to present the results for 
the (2, 1, s)-model. As the one-dimensional case with the 
theoretically assumed internal force projected onto the 
X-axis was already considered in the preceding subsec- 
tions, the analytical expressions for the first and second 
moments based on the two-dimensional Langevin equa- 
tion ((T|) are given by superposition of the motion in x- 
and y-direction. Thus, using Eqs. Q and ([5]) one obtains 
the vectorial mean position 



cos(0o) 
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(17) 



(r(i) - ro) = -(3FR' [l - e"^"-*] 
and the mean square displacement 

+ 2 (^/3^^') [Drt - 1 + e-^"*] . 

(18) 

As expected, the 0o-dependence vanishes in Eq. ([T51) due 
to the free translational motion in the two-dimensional 
plane. Furthermore, the diffusive term given by the first 
summand in Eq. (|18p is naturally twice as big as in 
Eq. ©. 



IV. ELLIPSOIDAL PARTICLE WITH ONE 
ORIENTATIONAL DEGREE OF FREEDOM 

We now generalize the previous considerations to el- 
lipsoidal particles. To cover the situation depicted in 
Fig. HKb) we have to take into account that, as opposed 



r 



to the case of spherical particles, the translational dif- 
fusion coefficient is anisotropic, which means that the 
diffusion tensor 



Dt = Daiii (E>u)+ Db{l - u «) u) 



(19) 



has to be applied. Here, I is the 2x2 unit matrix, 
u — (cos 4>, sin (j>) is the orientation vector, (g) a dyadic 
product, and Da and Di,, respectively, indicate the dif- 
fusion coefficients for translation in the direction of the 
two semi-axes of the ellipsoid. The index a stands for 
the semi-major axis, while b marks the semi-minor axis. 
Using the diffusion tensor (Eq. p^ ). the Langevin equa- 
tion for the center of mass position of the particle can be 
written in the form 



dr 

— = /3Dt • [Fu - VC/] + w . 



(20) 



Due to the anisotropy of the diffusion coefficient we can- 
not include the Gaussian white noise random force ex- 
actly in the same way as in Eq. (HI). Instead of that, 
we use the zero mean random noise source w(i). The 
variances of the components i,j € {x^y} are given by 
{wi{t)wj{t')) = 2Df{(j){t))6{t-t'). Thus, w,(t) are Gaus- 
sian random variables at fixed 4>{t). This follows the pro- 
cedure presented in Ref. 37] for "passive" ellipsoidal par- 
ticles. 



A. The (1, 1, e)-model 

It can easily be seen from Eq. (IT^ that Dt • u = 
DaU- Therefore, in the context of the (1, l,e)-model the 
Langevin equation for the center of mass position x of the 
self-propelled particle without external potentials can be 
written as 



Q.X 

-— = (3DaF C0S{(J)) +Wx 

at 



(21) 



The Langevin equation for the angle 0, which is needed 
in addition to Eq. I^T^ , does not differ from Eq. Q . For 
the following calculations, it is convenient to write the 



diffusion tensor (Eq. (IT91) ') as 
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where D = l/2{Da + -D&) marks the mean diffusion coef- 
ficient and AD the difference Da — Di, between the dif- 
fusion coefficients along the long and short axes of the 
ellipsoid. 

The mean position of an ellipsoidal particle is given by 
the first component of Eq. (|24l) . For the mean square 
displacement one obtains the relation 
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For an isotropic particle with AD = Eq. ([231) reduces 
to the result for a spherical particle in one dimension 
(Eq. ([5])), as expected. For an anisotropic particle the ad- 
ditional (second) term in Eq. (1^ yields a (/iQ-dependence 
at very early times, in the regime of bare translational 
diffusion, which is not present in the case of an isotropic 
particle (see the discussion in Sec. lIII At . This additional 
term represents the relative orientation of the initial di- 
rection of the long axis of the ellipsoidal particle and the 
direction of the linear channel. The (/)o- dependence can 
also be seen in Fig. [3 where the strength of the driving 
force is determined by the parameter F* = (3Fo \JDajDr 
for ellipsoidal particles. 

The analytical results for the third and fourth mo- 
ments, which are needed to calculate skewness and kur- 
tosis, are more complicated and, therefore, given in the 
appendix in Eqs. (|A.1[) and (|A.2[) . respectively. 



B. The (2, l,e)-model 

Again, we provide the analytical results for the case 
with two-dimensional translation as well. For an ellip- 
soidal particle, one obtains the expressions 
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for the first moment and 

{{r{t)-r,)')=4Dt + 2(^/3F^^ 

(25) 

for the mean square displacement, respectively. As in 
the mean square displacement of a spherical particle 
(Eq. (O), the (/)o-dependence also vanishes in Eq. (P5|. 
Furthermore, the contribution to the diffusive motion due 
to the initial orientation of the particle disappears for 
two-dimensional translation so that the diffusive motion 
is simply reflected by the term ADt in Eq. (I^5|) . 




FIG. 7. (Color online) Mean square displacement of an ellip- 
soidal particle with one orientational degree of freedom and 
anisotropy AD/Da ~ 0.8. Plots are shown for a) a self- 
propelled particle with F* — fiFoyjDa/Dr = 100 and b) a 
"passive" particle with vanishing effective force. The solid 
lines refer to a parallel and the dashed lines to a perpendicu- 
lar initial orientation. 



FREELY ROTATING SPHERICAL 
PARTICLE 



Up to now, particles with only one orientational de- 
gree of freedom have been examined. In this section 
we transfer our model to particles whose orientation is 
freely diffusing on the unit sphere. Considering a spher- 
ical particle this situation is shown in Fig. [Ijc). In 
the Cartesian lab frame the particle orientation u = 
(sin 9 cos (/?, sin 9 sin Lp, cos 9) is now given in terms of the 
two orientation angles 9 and (/?. Using the updated ori- 
entation vector the Langevin equation for the center of 
mass position is identical to Eq. ([1]) if the third compo- 
nent of all vectorial quantities is considered additionally. 

The orientational probability distribution for the freely 
diffusing orientation vector ^45.] is given by 

oo I 

P{9,^,t) = J2 E e-''^'^'+''>'Yr{9o,^o)Yr{0,^), 

1=0 m=-l 

(26) 
where YJ™ are the spherical harmonics. In Eq. ((^S)) we 
use the notation ^o = &{t = 0) and (po = f{t = 0) while 
the star indicates complex conjugation. 



A. The (1,2, s)-model 

To eliminate the (/3-dependence in the equation of mo- 
tion, we choose the z-axis to point in the direction of 
the linear channel, which we consider first. Moreover, we 
omit external potentials in the following. By calculating 
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FIG. 8. (Color online) Comparison of the mean position of a 
spherical particle with one and with two orientational degrees 
of freedom. The graphs for which the initial angle (po is given 
refer to the (1, 1, s)-model while the graphs designated by a 
certain value for the angle 9o show the results for the (1, 2, s)- 
model. In all cases the effective force is F* = 10 and the 
motion in z-direction is considered. 



{cos{9)) via Eq. (P5|) the first moment is obtained as 



{zit)-za) = -PFR^cos{0o){l 



-2Dr-t\ 



(27) 



The mean position in the (1, 2, s)-model is very similar 
to the same in the (1, 1, s)-model (see Eq. 2]). Here, the 
azimuthal angle 9 takes the role of the angle cj) in the 
(1,1, s)-model. In particular, the two results agree up to 
linear order in time i, whereas they deviate for longer 
times due to the enhanced probability of the sphere with 
full orientational freedom to assume a configuration with 
an orientation pointing in the direction of the equator. 
This, in turn, on average causes a smaller force com- 
ponent along the z-axis and a smaller plateau-value of 
the excursion limt_j.oo(^(^) ^ ^(0)), which is illustrated in 
Fig.H 



Using Eq. ((26|) and, thus, the fact that every function 
that depends exclusively on 9 and ip can be expanded as a 
linear combination of spherical harmonics, we also obtain 
the analytical result for the mean square displacement. 



which is given by 
{{z{t)-zor)^lR'Drt+(^PFR' 
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(28) 



Comparing Eqs. (|28l) and (O (as shown in Fig. ^ it 
turns out that also the mean square displacements of the 
spheres with two or one orientational degree of freedom 
in a linear channel are almost identical. In particular. 
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FIG. 9. (Color online) Comparison of the mean square dis- 
placement of a spherical particle with one and with two ori- 
entational degrees of freedom. The remarks in the caption of 
Fig. [8] are also valid for this figure. 



their functional forms are the same up to second or third 
order in time t, for the cases of a parallel (0o = ^o ~ 0) 
or a perpendicular (0o = ^o ~ '^/'^) initial configuration, 
respectively. Therefore, the crossover time scale from 
the diffusive to the super-diffusive regime ti is exactly 
the same as in the (1, 1, s)-model and given by Eq. ()10|) . 
The second time scale from the super-diffusive to the 
diffusive regime is in the (1, 2, s)-model given by ^2 = 
{2Dr)^^ and therefore half as large as in the (l,l,s)- 
model. Concomitantly, the long-time diffusion constant 
is smaller as compared to Eq. (jlip and given by 



Dr = -DrR^ 



l + -i(3RFy 



(29) 



As for the long-time limits of the mean positions, the 
difference to Eq. ([TH) reflects the fact that the freely ori- 
ented sphere is more likely oriented perpendicular to the 
channel axis than the sphere that is confined to rotate in 
the plane. 

The non-Gaussian behavior of the self-propelled parti- 
cle with free orientation on the unit sphere is embodied 
in the third moment 
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and in the fourth moment 
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The curves for the skewness (see Fig. [TU|) and the kurto- 
sis (see Fig. [TT|) of the probabihty distribution function 
^{x,t) for the (1, 2, s)-model are obtained by shrinking 
their counterparts for the (1, 1, s)-model in x-direction as 
well as in the direction of the t-axis. This is very similar 
to the findings concerning the mean position of the par- 
ticle (see Fig. [8]) . The extrema of skewness and kurtosis 
and the change of sign of the kurtosis, that is observed for 
non-perpendicular initial configurations, already occur at 
smaller times. Obviously, the existence of the negative 
1/t long-time tail in ^{t) (see Sec. IIII A[) is not affected 
by the number of orientational degrees of freedom of the 
particle. 



B. The (2,2, s)-model 

For more than one translational degree of freedom, if 
the particle motion takes place in the xy-plane, the first 
and second moments are given by 
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FIG. 10. (Color online) Comparison of the skewness 5(4) 
of the probability distribution function '^(x,i) for a spheri- 
cal particle with one and with two orientational degrees of 
freedom. The solid and the dash-dotted lines refer to the 
(1, 1, s)-model while their counterparts for the (1, 2, s)-model 
are given by the dashed and the dotted lines. 



respectively. As before, the results for the mean position 
and the mean square displacement are almost identical 
with respect to their lower-dimensional counterparts. Be- 
side the change from a cosine to a sine in the last part of 
Eq. ([55]) . which is due to the change of accessible dimen- 
sions, the only difference of the mean square displace- 
ments in the (2,2,s)- and in the (1, 2, s)-model consists 
in an additional factor of 2 for all terms that do not de- 
pend on ^0- 
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FIG. 11. (Color online) Comparison of the kurtosis 7(f) of 
'i'{x,t) for a spherical particle with one and with two orienta- 
tional degrees of freedom. See also the remarks in the caption 
of Fig. \W\ 



C. The (3, 2, s)-model 

The (3, 2, s)-model is the most general situation con- 
cerning a spherical particle. In this case, the mean 
position of the particle is obtained by adding the z- 
component (Eq. (P7)) ) to Eq. (P^ . The mean square 
displacement is given by 



((r(t) - ro)2) = SR^Drt 



\ {>-' 



[2Drt - 1 + e-^-D'-*] 



(34) 



The simplicity of Eq. (|34p can be explained by the fact 
that no dependence on the initial orientation can appear 
due to the completely free motion |53j . 



VI. FREELY ROTATING ELLIPSOIDAL 
PARTICLE 



To complete our examination of the different model sit- 
uations, we now consider a freely rotating self-propelled 
ellipsoidal particle as sketched in Fig.[ljd). The (£>, 2, e)- 
models, where D E {1,2,3} is the number of transla- 
tional degrees of freedom, are based on the Langevin 
equation (|20p . To transfer this two-dimensional equation 
to the three-dimensional case regarded here, the orienta- 
tion vector u = (sin 9 cos (p, sin 9 sin (f, cos 9) is used and 
the third component is added to all vectorial quantities. 
The explicit form of the diffusion tensor Dt is given by 
Eq. (J19p by means of the updated orientation vector. 



A. The (1,2, e)-model 

As in the previous sections, we begin by considering 
the case of one-dimensional translational motion with- 
out external potentials, i.e., with the (1, 2, e)-model. Us- 
ing the ansatz based on spherical harmonics according to 
Eq. (P5|) for an ellipsoidal particle leads to the analytical 
results 



(z(i) - zo) = \l3F^ [I - 6-2^--*] cos(0o) (35) 



for the mean position and 
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(36) 



for the mean square displacement, respectively. Corre- 
sponding expressions for the third and fourth moments 
and, thus, for skewness and kurtosis were calculated as 
well. As these are quite lengthy, they are given in the ap- 
pendix in Eqs. (IA.3I) and (IA.4I) . respectively, but are pre- 
sented graphically in Figs.[T2]and[T3l Different regimes of 
non-Gaussian behavior are manifested by different signs 
of the kurtosis. For "passive" particles with vanishing 
internal effecitve force F* — (solid line in Fig. [T5)) no 
change of sign is observed. The kurtosis is positive and 
a simple maximum occurs at t ^ t2 = {2Dr)~^. These 
findings correspond to the results for "passive" ellipsoidal 
particles in two dimensions that were studied in Ref . [33] . 
In contrast to simple non-Gaussian behaviour, the situa- 
tion turns out to be much more complex if self-propelled 
particles are considered. If the self-propulsion outweighs 
the effect of the kicks of the solvent particles (dashed line 
and dotted line in Fig. [T51) , several maxima, minima, and 
changes of sign induce a rich structure in 7(t) indicating 
different non-Gaussian behavior at different time scales. 
The characteristic 1/i long-time tail observed for spher- 
ical particles is also found for ellipsoidal self-propelled 
particles. The skewness S{t) (see Fig. [T2|), which is a 
measure of the asymmetry of the probability distribu- 
tion, reveals a higher degree of complexity as well. While 
S{t) is zero for "passive" particles (solid line in Fig. [T2|) . 
this parameter also shows a much richer structure if self- 
propelled particles are regarded. 
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FIG. 12. (Color online) Skewness S{t) of the probability dis- 
tribution function '^{x,t) for an ellipsoidal self-propelled par- 
ticle as a function of time. This figure illustrates the analytical 
results of the (1, 2, e)-model. While the initial angle Oq and 
the anisotropy ^D jDa are constant as given in the figure, 
graphs for various values of the effective force Fl are shown. 
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FIG. 13. (Color online) Kurtosis 7(t) of *(a;,t) for an ellip- 
soidal self-propelled particle as a function of time. The graphs 
correspond to the graphs in Fig. [T2]as far as the values of the 
various parameters are concerned. 



B. The (2, 2,e)-modeI 



By simply combining the results for one-dimensional 
translation in x- and in j/-direction, for the (2, 2, e)-model 
we obtain the expressions 
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for the mean position and 
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for the mean square displacement. 



C. The (3, 2,e)-model 

In a last step, we now add the third translational de- 
gree of freedom. Hence, in this subsection we consider 
the most general case with free translational and rota- 
tional motion of an ellipsoidal self-propelled particle. By 
adding the third component (Eq. ([55])) to Eq. ([57]) we 
obtain the result for the mean position. The analytical 
expression for the mean square displacement is simply 
given by 

((r(i) - ro)^) = &Dt - ADt 



l(,Ff ) [2D,,-1 



-2Drt-\ 



(39) 



(37) 



To explain the simplicity of this result we want to point 
to the short discussion after Eq. ([M)) . 



VII. CONCLUSION 

In conclusion, we have analytically solved the Brown- 
ian dynamics of an anisotropic self-propelled particle in 
different geometries by presenting explicit results for the 
first four moments of the probability distribution func- 
tion for displacements. The particle is driven along an 
axis which itself fluctuates according to rotational Brow- 
nian dynamics. After a transient regime which is char- 
acterized by two distinct time scales, there is diffusive 
behavior for long times. The results for the long-time dif- 
fusion constants D^ for the different groups of model sit- 
uations are given in TableUl For intermediate times, non- 
Gaussian behavior is revealed by a non- vanishing kurtosis 
in the particle displacement which decays as 1/t for long 
times t. For special initial conditions (nearly perpendic- 
ular initial orientation and large effective forces) , we find 
a superdiffusive transient regime where the mean square 
displacement scales with an exponent 3 in time. The ana- 
lytical results can be used to compare with experimental 



12 



TABLE I. Long-time diffusion constant Dl for the different 
model situations. Using the definition in Eq. (|lip the analyt- 
ical results are directly obtained from the respective results 
for the mean square displacement. 



Model 



Long-time diffusion constant 



{D, 1, s) _ {i/3)DrR^ [1 + {2/3){pFRf] 

{D, 1, e) D + l/{2Dr) {PF [D + {l/2)AD)f 

{D, 2, s) _ (4/3)D,i?2 [1 -f (2/9)(_^Fi?)2] 

(D, 2, e) D~ (1/6) AD + l/{&Dr) (/3F (D + (1/2) AD) )^ 



systems of, e.g., swimming bacteria or self-propelled col- 
loidal particles. Any deviations point to the importance 
of hydrodynamic interactions with the substrate and with 
neighbouring particles at finite density. 

It would be interesting to generalize the analysis to- 
wards various situations. First of all, hydrodynamic in- 
teractions were neglected in our studies. While this is 
justified in the bulk, hydrodynamic interactions become 
important at finite densities [5J| and may significantly 
influence the distribution of the mean square displace- 
ments 



55| . If the particle is moving close to a substrate 



and the substrate is hydro dynamically not flat, then hy 
drodynamic interactions play a significant role, too [56{ . 
Second, it would be interesting to include an additional 

I 



torque in the Langevin equations of motion. This leads 
to circular motion in two dimensions [32l | while for three 
spatial dimensions helical motion is expected as also sug- 
gested by a slightly different model presented in Ref. fS^ . 
Next, a non-Gaussian noise [58] in the Langevin equa- 
tions might be relevant for modeling real swimming ob- 
jects [13|,|57[. Furthermore, the actual propulsion mecha- 
nism was modelled just by an effective force. A consider- 
ation of more details about the actual propulsion mecha- 
nism might be necessary to analyse short-time dynamics 
in more depth. In addition to that, the model migth be 
transferred to more complicated geometries [59* such as 
ratchets |^6(i 6l|, for example. Finally, the collective 



dynamics |62h64| of many swimmers will lead to further 
effects like swarming, swirling and jamming. While at 
high densities hydrodynamic interactions are expected 
to play a minor role, the direct particle-particle interac- 
tions get relevant and should be incorporated in theory 
[65l | and simulation ;66i] . 
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Appendix: Further analytical results 



In the following, we summarize the analytical results for the third and fourth moments of an ellipsoidal particle. 
For reasons of clarity, we use the notation F* = jiFQ^Da/Dj.^ the ratio 5 = /S.D/Da and the scaled time r = Drt. 
For the (1, 1, e)-model the third moment is given by 



{x(t) - Xpf 



F* 



-— cos(0) -f — cos(3(?!)) -t- — cos((/))e ^ — — cos(3 0)e "^ -I- 3 cos(0)t — — cos(0)e ^'^ 



40 



24 



cos(3(/))e 



cos((/))Te 



240 



16 



e-^^cos(3 0) 



6^; (1--5 



cos(0)r — cos(0)r e "^ 



3F:6 



5 2 _ 

cos{(f>) + — cos{3(j)) - - cos(0)e ^ 



16 



cos{3(j))e "^ 



13 

720 



-9t 



1 



cos(3 '/') - 2 cos(<?!))t e ^ + ^ cos(0) 



-At 



24 



40 



cos{3 (j))t 



-At 



(A.l) 



and the fourth moment is 
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The corresponding results for a freely rotating elhpsoidal particle ((1,2, e)-model) are the third moment 
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+ ^ cos(0)e 



-12 r 



1 

90 



1 



e-'''^cos(3 6') eos(3 

^ ^ 40 



-2t 



1 

10 



cos(6')r e 



-2r 



(A.3) 



and the fourth moment 
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l/3Ff + AFf (1-1/25) -4S {I- 1/2S) 

295147905179352817853 , ^, 

442721857769029238784 '^ . 

'pp*2c/ 16971004547812787156293 184467440737095508157 _q^ 
" V 79689934398425262981120 ^ 55782954078897684086784*^ 



184467440737095508157 



18594318026299228028928 
230584300921369395 



cos(2 6')e 



-6t 



9223372036854775808 



cos(2 6')e" 



230584300921369395 
' 9223372036854775808 
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-2t 
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cos(2 6l) 



i6' 



-l/9cos(2.) + ^-^e--leos(2.)e- 



F- ( I eos(2 .)e- + 1/6 eos(2 .) - I e- - ^ e- - 1 cos(2 .)e- - g 



-6t 



+ 1/6 i^r (1 - 1/2 (5) -10 + 6 cos(2 6i) + 9 e"^ ^ - 9 cos(2 6i)e-^ ^ + 3 cos(2 e)e 



-6t 



+ 2/3 (5(1- 1/2 (5) 3 cos(2 6*) + 1 - 3 cos(2 6i)e 



-6 r ^^6 r 



^* 



2400 



cos(2 6')e 



-12t 



-12t 



9600 



672 



cos(4 6')e 



-6t 



441 



cos(2 6')e 



'6t 



— cos(2 61) - — cos(4 e)e-^ ^ + — cos(2 6')e-2 ^ + ^ cos(4 g) 
720 ^ ^ 480 ^ ^ 300 ^ ^ 960 ^ ' 



11 



4901 



1411740617885935 



14112 

1 
23520 ' 



14400 73786976294838206464 



273 
800' 



-20 T 



cos(2( 



1 



1920 



cos(4 6')e 



-12t 



1 



-20 T 



cos(4 6i) 



QFrS 



304002342334733405771 

2361183241434822606848 ^ 

26235369349275805845 

" 18889465931478580854784 
491913175298921373631 
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298099384231 14633542725 



13440 
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318759737593701051924480 



cos(2 61) 



3123845428418270308859904 
271536072765004556323 



cos(4 6')e" 
cos(4 61) 
cos(2 6l)e" 



2^ 5247073869855161169 _i2r 
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24595658764946071489 ,,^, .^ 

cos(4 6')e ^^ 
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807229520665529965278287 
" 6375194751874021038489600 



21250649172913403461632 
2951479051793528185637 



cos(2 6')e" 



807045053224792883 

' 166020696663385964544 
2951479051793528185637 



cos(4 6i)e"2^ 

„-20t . 



23139595766061261547110400 
5247073869855161169 
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cos(2 6l)e 
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'cos(2( 



"eos(4 0) 
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Z5' 



-20 T 
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cos(4 6') 



-20 r 
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008(2 6*) 



84 



cos(4 6')e 



-6r 



19 



37 



1800 1323 



cos(2 6')e 



-6t 



7 



17 



eos(4 0) + ^ eos(2 0) 

120 ^ ^ 270 ^ ' 1764 



-6r 



-20t 



9800 



(A.4) 
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